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The frequencies of transverse-optical V phonons in KNDO3 and KTa03 are calculated in the 
frozen-phonon scheme making use of the full-potential linearized muffin-tin orbital method. The 
calculated frequencies in the cubic phase of KNbC>3 and in the tetragonal ferroelectric phase are in 
good agreement with experimental data. For KTa03, the effect of lattice volume was found to be 
substantial on the frequency of the soft mode, but rather small on the relative displacement patterns 
' ■ of atoms in all three modes of the T\ u symmetry. The TO frequencies in KTaC>3 are found to be of 

\ the order of, but somehow higher than, the corresponding frequencies in cubic KNb03. 
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I. INTRODUCTION 

\o : 

Evaluation of vibrational properties within the frozen-phonon scheme is known to be a hard test on the quality 
of full-potential total-energy calculations. The curvature of the total-energy surface over the manifold of various 
atomic displacements is much more sensitive to the details of the calculation scheme than merely the position of the 
• total-energy minimum, i.e., the equilibrium geometry. However, in case of success (that can be easily checked by 
comparison with experimentally measurable phonon frequencies) the calculation provides substantial information on 
the microscopical driving forces behind the specific vibration patterns and may give some insight into the dynamic 
properties of the crystal in question. For ferroelectric materials, the interest in the phonon calculations is motivated 
by the apparently crucial role of zone-center phonon softening in the mechanism of ferroelectric phase transition which 
is subject to controversial discussions 0J|,^,Q. 

An early calculation of phonon dispersion curves within the empirical shell model has been done for KNbC>3 
by Fontana et al. jjj. First-principles calculation of V phonons in BaTiC>3 which has the same crystal structure 
and exhibits the same sequence of ferroelectric transitions as KNbC>3 have been done by Cohen and Krakauer [||. 
Liechtenstein et al. |?J analyzed two vibrational modes (however not related to ferroelectric transition) in another 
perovskite - BaBiC^. Singh and Boyer calculated V and R phonons in cubic KNbC>3. Recently, Zhong et al. |J 
£h , obtained TO and LO T-phonon frequencies in a number of cubic perovskite-type ferroelectrics. 

In the present paper, we continue the comparative research of KNbC>3 and KTa03 initiated in the previous ab initio 
study of equilibrium geometry (Ref. pc| ], referred to further as I). In the present paper, we concentrate on phonon 
^ • frequencies in these compounds making use of the same calculation scheme (full-potential linear muffin-tin orbital 



code by Methfessel [11 LJ) and setup as discussed in I. 

For KNb03, we performed the T-phonon calculations for the nonpolar cubic phase and for the first (as the temper- 
ature lowers) ferroelectric phase, i.e., the tetragonal one. The results are discussed in Sees. || and III. In Sec. IV 
the results for KTaC>3 in the cubic phase for two lattice spacings are presented. 



II. KNbO s : CUBIC STRUCTURE 



As has been pointed out in I, the calculated (from the total-energy minimum) equilibrium volume of the cubic 
phase of KNb03 turns out to be ^95% of the experimental cell volume (extrapolated to zero temperature). Such 
discrepancy is known to be typical for calculations based on the local density approximation (LDA). The effect of 
different volumes on the trends in total energy lowering and related equilibrium displacements is discussed at length 
in I. Since the curvature of the total-energy hypersurface may be affected by the error of about 5% in the cell volume, 
we preferred to perform our phonon calculations for the experimental lattice spacing, in order to produce results 
better comparable with experimentally measured phonon frequencies. Specifically, we took a=4.00 A for the cubic 
phase (extrapolated value on T=0 from lattice constants of the high-temperature phase, according to Ref. |ll|). The 
measured lattice constant in the cubic phase which exists above 418 °C is about 4.02 A, slightly increasing with 
temperature. [jl3f The variation of volume over the temperature range is therefore much smaller than the difference 
between experimental and theoretical volumes and is ignored in the present study. 
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As it is known (see, e.g., Refs. [pi and [H), T-phonon vibration modes in the cubic perovskite structure, after 
projecting out three translational modes, are split by symmetry into three T\ u modes and one T 2u mode (all of 
them are triple degenerate). We followed essentially the guidelines of Cohen and Krakauer || in projecting out the 
translational modes, as well as in the choice of appropriate symmetry coordinates. When constructing the force 
constant matrix, we performed about 30 calculations for different displacement patterns, in order to provide a good 
multidimensional fit for the total-energy hypersurface by a fourth-order polynomial, retaining then appropriate second 
derivatives. In such a way, we could obtain stable and controllable values, e.g., for off-diagonal elements of the force 
constant matrix. 

The simplified procedure of extracting each of these elements independently from a single calculation involving 
particular combined displacement has been used by Cohen and Krakauer Q and applied to KNb03 by Singh and 
Boyer Q . In our opinion, this scheme may be not sufficiently reliable in case of an essentially nonquadratic shape of 
the total-energy hypersurface, resulting in the apparent dependence of the calculated force constants on the particular 
displacement chosen. 

On solving the secular equation of the 3x3 lattice dynamics problem (for the T\ u mode), 

[GF - ll> 2 }u = 

(Rcfs. |l5| and |l6|), with G and F matrices defined as in Ref. Q, the pattern of Cartesian displacements of individual 
atoms is restored by back symmetry transformation. For the atoms with the fractional coordinates in the cubic 
perovskite cell K(0, 0, 0); Nb(0.5, 0.5, 0.5); Oi(0, 0.5, 0.5); O 2 (0.5, 0, 0.5); O 3 (0.5, 0.5, 0), the vibrations, e.g., along 
the [001] direction are given by 

/zi\ / 4 -1 

z 2 -14 

z 3 = -1 -1 

z 4 -1 -1 

\zj \ -1 -1 

These displacements related to the center of mass and multiplied by square roots of individual atomic masses produce 
the orthogonal eigenvectors of vibrational modes, which are presented, along with the frequencies, in Table | in 
comparison with some experimental data. 

Our calculations reproduce the measured frequencies of TO2 - TO4 phonon modes reasonably well, while the TOi 
mode frequency is imaginary, as can be expected for a soft mode at zero temperature. In accordance with our previous 
analysis of I and in agreement with the results of Singh and Boyer, |^[ the eigenvector corresponding to the soft mode 
represents roughly a displacement of Nb with respect to the rest of the crystal. It should be noted that in I, the 
off-center displacement of Nb was found to be energetically favourable not only at the experimental cell volume for 
a =7.553 a.u., but also at the somewhat smaller optimized theoretical volume (a =7.425 a.u.). Therefore, the soft 
mode frequency, even if calculated at the latter volume, is expected to be imaginary. On the contrary, Singh and 
Boyer || could only obtain the energy lowering by Nb displacements at a =7.589 a.u. (the value taken in Ref. ||] 
for the experimental lattice constant), but not at their theoretical lattice constant of 7.448 a.u. This discrepancy is 
probably related to some differences in the calculation schemes and needs to be further investigated. 

Our results indicate a good agreement with the experimental data in determining the frequency of the TO2 mode, 
whereas the frequencies of TO3 and TO4 modes are systematically underestimated in both the present calculation 
and that of Singh and Boyer. Q The reason may be that the TO2 mode is essentially the vibration of K atoms 
against the rest of crystal. As soon as the potential well related to individual off-center displacements of K is the 
most parabolic comparing to potentials felt by all other atoms (it has been discussed at length in I), the harmonic 
approximation seems to work best for this kind of displacement. Additional evidence confirming this point is that the 
experimental frequencies measured for this mode are most stable over a broad temperature range. The reason why the 
corresponding frequency is underestimated by 15% in the calculation of Singh and Boyer ||] may be the unadequacy 
of their fit for the total-energy hypersurface which has been constructed only on six different displacement patterns 
to span all three Ti^-type modes. 

TO3 and TO4 modes involve two different kinds of stretching of the oxygen sublattice which seem to give rise to a 
rather unparabolic total-energy surface, so the accuracy of the phonon description in the harmonic approximation is 
not sufficient. 

Turning to a more detailed analysis of the eigenvectors, one may conclude that there is an overall agreement between 
the displacement patterns calculated in the present work and in Ref. ||, for the experimental lattice constant. (There 
is also clear similarity with the displacement patterns calculated for T\ u modes of BaTiOa by Cohen and Krakauer ||). 
However, in our calculation the displacement of K atoms in the soft mode and in the TO3 mode is more pronounced. 
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The displacement pattern within the soft mode, although calculated from the second derivatives of the total- 
energy hypersurface in the ideal cubic structure, contains essential information about relative finite displacements of 
atoms in the course of cubic to tetragonal phase transition. This has been shown already by Cohen and Krakauer 
H for BaTi03. For KNbC>3, the relative displacements corresponding to the eigenvector of the soft mode in our 
calculation are presented in Fig. ^(a) in comparison with the atomic positions in the tetragonal cell as determined by 
neutron diffraction jl8| [Fig.[l](b)]. Whereas the experimental coordinates may be interpreted as being due to opposite 
movement of Nb and O sublattices, with the K atoms relatively undisplaced from their positions in the center-of-mass 
scale, our calculations show a pronounced tendency of K atoms to remain stuck to the oxygen sublattice, so that 
the ferroelectric transition may be roughly considered as primarily due to net off-center Nb displacement versus the 
rest of the crystal. It should be noted that the K displacement has been experimentally estimated with the lowest 
accuracy, compared to other atoms, as has been pointed out in Ref. |18|. 



III. KNbOa: TETRAGONAL STRUCTURE 



For the tetragonal phase, the situation is complicated by the presence of strain. In I, we have shown that the 
tetragonal strain can be optimized, along with the off-center displacement of Nb atoms and under the constraint 
of constant (theoretical) cell volume, to be in fairly good agreement with the experimental estimate of c/a ^1.02. 
Now that we would like to proceed with the experimental lattice constant which does not provide the minimum of 
the calculated total energy, the optimization of the strain from first principles is no longer justified. Instead, we 
should take the experimental lattice constants at some temperature as an external constraint. The choice a=4.00 A, 
c/a=1.0165 corresponds to what is measured at about 250°, near the lowest-temperature range of existence of the 
tetragonal phase. [13| Within thus predetermined crystal lattice, associated with any one type of atom (for instance, 
Nb), K and two nonequivalent types of O are free to relax along the [001] axis towards equilibrium positions compatible 
with tetragonal symmetry. Since this equilibrium geometry is not determined by symmetry as in the cubic phase, it 
must be found from first-principles calculations prior to further frozen-phonon analysis. This search for the global 
total-energy minimum over three independent parameters has been accomplished by a polynomial fitting and needed 
several tens of total-energy calculations to achieve sufficient accuracy. 

The optimized perturbations of the fractional coordinates of atoms in the cubic perovskite cell, accounting for 
the off-center displacements along [001] in the tetragonal structure, are found to be the following: 0.046(K); 0(Nb); 
0.034(Oi, O2); 0.047(03). The arbitrary choice of unshifted Nb sublattice is taken here in order to enable a direct 
comparison with the experimental data of Hewat |l^] which are correspondingly 0.018; 0; 0.040; and 0.044. The 
absolute displacements of atoms from their symmetry positions are shown in Fig. |l|(c) , where the center of mass has 
been kept fixed. 

As was the case with the soft mode analysis in Sec. Q, K shows, according to our calculation, a pronounced 
tendency to remain stuck to the oxygen sublattice. This is identical to what we have found in I for the rhombohedral 
phase, when considering the energetics of coupled K and O distortions. Otherwise, the main trend in the ferroelectric 
transition, namely that the largest displacement is that of Nb with respect to O3, as well as the magnitude of this 
displacement, are in fairly good agreement with the experiment. 

As the crystal space group is reduced from PmZm to PAmm at the ferroelectric phase transition, the vibrations 
along the [001] direction are no more degenerate; three corresponding T\ u modes which retain the tetragonal symmetry 
of the crystal lattice belong now to the A\ representation, whereas the formerly T^u mode now becomes B\. After 
projecting out the uniform displacement along [001], we arrive at the following symmetry coordinates (similar to those 
of Ref. §): 
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Vibrations normal to the tetragonal axis remain doubly degenerate; former T± u and T2 U modes however are now 
mixed in the E representation, with the following symmetry coordinates (along [100]): 
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and analogously along [010]. 

Calculated phonon frequencies and restored orthogonal eigenvectors (individual atomic displacements, multiplied 
by square roots of masses) are presented in Table 0; a comparison with the experimentally determined frequencies is 



shown in Table III. With the exception of the soft TOi mode, experimental frequencies are fairly stable within the 
temperature range of the tetragonal phase, and calculated frequencies are in all cases in reasonable agreement with 
them. The best agreement is obtained for the modes involving a relatively small amount of Nb vibration, which is 
known from our previous analysis in I to be essentially anharmonic. The most parabolic potential well, as has been 
already emphasized in Sec. [n], is that related to the K atom; as a result, the harmonic frequency of the TO2--E mode, 
which represents roughly the vibrations of K against the rest of crystal, is in fairly good agreement with experiment. 
The same applies to the TO3 — A\ mode which is essentially the vibration of the basal O atom against all others. 
Consistent with this point of view, the experimental frequencies for these particular modes are the most stable over 
the temperature. Other modes either include a considerable contribution of the Nb displacement (two lowest A\ 
modes), or represent essentially the stretching of the oxygen cage (B\ and two highest E modes), in which case the 
harmonic approximation is less accurate, and the temperature dependence of frequencies should be noticeable. 

The double degenerate soft mode exhibits the tendency of atoms to shift from their positions, which were optimized 
with respect to displacements along [001], but correspond to a saddle point of more general total-energy hypersurface 
(see I) in one or another orthogonal direction, compatible with orthorombic crystal structure. The eigenvector of the 
soft mode indicates that such transformation involves primarily the displacement of Nb sublattice against the rest of 
crystal, as was the case in the cubic to tetragonal phase transition. 

IV. KTa0 3 

KTaOs, in contrast to KNb03, does not undergo a ferroelectric phase transition, remaining in the cubic perovskite 
structure over the whole temperature range. This crystal seems to be however at the very threshold to a ferroelectric 
phase transition, as is indicated by considerable softening of its TOi mode at low temperatures (see, e.g., Refs. ||l]] 
and [Q) and by the fact that this transition can be induced by applied uniaxial stress. [^3| In I, we studied the 
energetics of off-center displacements of different atoms in KTa03 and found the corresponding potential well to be 
the most anharmonic and volume-dependent for Ta. We found the cubic phase to be stable at the theoretical (i.e., 
underestimated by ^5% as is typical with the LDA) volume, but unstable towards off-center Ta displacements at 
the experimental volume. In order to study in more detail the effect of combined atomic displacements near the 
ferroelectric threshold, we calculated the phonon frequencies and eigenvectors at both theoretical (a=3.928 A) and 
experimental (a=3.983 A) values of the lattice constant. The calculation setup (basis set, choice of radii) was the 
same as described in I; the symmetry analysis for phonons is identical to that in the cubic phase of KNbC>3 . 



The calculated results are shown in Table IV in comparison with the frequencies measured by several techniques. 
Some more experimental data, which however fall within the same limits, may be found in Refs. pH and p7| . 
The above-mentioned tendency of Ta to go off center, as found in the calculations performed for the experimental 
lattice constant, results in imaginary frequency of the soft mode, the displacement pattern of which is essentially the 
movement of Ta against the rest of crystal. At the theoretical volume, which is smaller, the soft mode frequency was 
found to be real, numerically close to the data experimentally measured at about 200 K, and consequently higher than 
the experimental results obtained at lowest temperatures. The real behavior of the soft mode seems to fall within the 
limits provided by these two cases. 

In spite of the differences in the calculated soft mode frequency, the corresponding relative displacements of atoms 
determined for both lattice constants are very close. One can note however the difference in the vibration of the 
oxygen sublattice. At the theoretical lattice constant, i.e., for the case of a nonpolar phase and real frequency, the 
oxygen cage remains rigid on the vibrations; as the volume increases and the ferroelectric transition develops, O3 
atoms shift further towards Ta ions than O1.2 atoms do. The resulting slight distortion of the oxygen cage is typical 
for ferroelectric KNb03 according to our calculation (Table ||) and for BaTi03. || Similar trends are seen in the 
calculation by Singh and Boyer || for two lattice constants of KNb03 . 

The volume dependence of the frequencies of two higher T± u modes, although much less dramatic than it was for 
the soft mode, is nevertheless noticeable. The reason for the fact that frequencies at larger volume are systematically 
lower is that the potential wells related to off-center displacements of all constituent atoms become broader as volume 
increases, with smaller curvature at the equilibrium, as has been studied in I. 

We cannot provide an explanation for the opposite tendency found for the T^ u mode. It seems not to be an artifact 
of our calculation because Singh and Boyer Q obtained the same trends in their calculations for two lattice constants 
in KNb0 3 . 
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The whole set of calculated frequencies lies slightly but systematically higher than those in KNDO3. KTaC>3 therefore 
can be regarded to be the stiffer crystal as is also evident from its higher bulk modulus (calculated in I) . 

V. SUMMARY 

We calculated frequencies and eigenvectors of TO T phonons in cubic and tetragonal KNbC>3 and in cubic KTaC>3. 
For cubic KNbC>3, the soft T\ u mode is characterized by an atomic displacement pattern which is close to the pattern 
of the off-center displacements in the ferroelectric tetragonal phase, as determined by explicit structure optimization. 
The frequencies of the other two T\ u modes and the T 2u mode are in good agreement with experimental data. In 
the tetragonal phase, the softening of an E mode indicates an instability towards the transition to the orthorombic 
ferroelectric phase. The frequencies of the other six A 1; Bi, and E modes are found from the calculation to be in 
reasonable agreement with experiment - especially for modes whose atomic displacements mostly lie in quasi-harmonic 
regions of the energy hyperspace. In KTaC>3, the TO frequencies and eigenvectors have been calculated for two values 
of the lattice constant in the cubic phase. Although the tendency for the ferroelectric transition and correspondingly 
the calculated frequency of the soft mode are found to be quite sensitive to the volume, the displacement patterns 
within the soft mode and other T iu and T 2u modes arc not. The frequencies are found to be close to, but systematically 
slightly higher than, the corresponding frequencies in KNb03. 
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FIG. 1. Relative displacements of atoms in KNb03 along [001]: from the soft mode ei gen vector in the cubic structure [(a), 
arbitrary scale]; from neutron diffraction measurements in the tetragonal structure, Ref. (lq| (b); from geometry optimization 
in the tetragonal structure (c). 



TABLE I. Calculated T-TO frequencies and eigenvectors in cubic KNbO;?. 
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a for a=4.016 A. 

b Infrared reflectivity measurements at 710 K, Ref. [14|. 
c Infrared reflectivity measurements at 740 K, Ref. [ 1 7 L 
d Infrared reflectivity measurements at 1180 K, Ref. |14|]. 

e Infrared reflectivity measurements at 585 K (in the tetragonal phase), Ref. jF 



TABLE II. Calculated F-TO frequencies and eigenvectors in tetragonal KNbOg. 
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TABLE III. Calculated and measured frequencies of F-TO phonons in tetragonal KNb03 
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a Infrared spectroscopy at 585 K; Ref. fll - 

b Neutron spectroscopy at 518 K (Ref li*9l0 ; numerical values cited in Ref. [|20|. 
c Raman spectroscopy at 543 K; Ref. |2C 
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TABLE IV. TO frequencies and eigenvectors in KTaC>3. 
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a Infrared reflectivity measurements at 12—463 K, Ref. ]24j ]. 

b Hyper- Raman scattering measurements at room temperature, Ref. p5|. 

c Raman scattering measurements at room temperature (soft mode) and at 10 K, Ref. p 
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